Possible evidence of extended objects 
inside the proton 



R. Petronzio 1 , S. Simula 2 and G. Ricco 3 

x Dip. di Fisica, Universita di Roma "Tor Vergata" and INFN, Sezione di Roma II, 
Via della Ricerca Scientifica 1, 1-00133 Roma, Italy 
2 Istituto Nazionale di Fisica Nucleare, Sezione di Roma III, 
Via della Vasca Navale 84, 1-00146, Roma, Italy 
3 Dip. di Fisica, Universita di Genova and INFN, Sezione di Genova, 
Via Dodecanneso 33, 1-16146, Genova, Italy 



Abstract 

Recent experimental determinations of the Nachtmann moments of the inelastic struc- 
ture function of the proton F%(x, Q 2 ), obtained at Jefferson Lab, are analyzed for values 
of the squared four-momentum transfer Q 2 ranging from ~ 0.1 to ~ 2 (GeV/c) 2 . It is 
shown that such inelastic proton data exhibit a new type of scaling behavior and that 
the resulting scaling function can be interpreted as a constituent form factor consistent 
with the elastic nucleon data. These findings suggest that at low momentum transfer 
the inclusive proton structure function originates mainly from the elastic coupling with 
extended objects inside the proton. We obtain a constituent size of « 0.2 -j- 0.3 fm. 



PACS numbers: 13.60.Hb, 14.20.Dh, 13.40. Gp, 12.39.Ki 



1 Introduction 



Since long time hadronic spectroscopy and Deep Inelastic Scattering {DIS) data have been 
the two main sources of information on hadron structure. The investigation of hadron mass 
spectra has led to the introduction of the concept of quarks [T|, leading to the very fruitful 
idea that meson and baryons are bound-states of two and three quarks. Such quarks are 
commonly referred to as Constituent Quarks (CQ's). The DIS data (starting from the 
pioneering experiments at SLAC in the sixties j2]) have been successfully interpreted in 
terms of a short-distance partonic structure of the hadrons, i.e. the presence of point-like 
constituents inside the hadrons 

With the advent of Quantum ChromoDynamics (QCD) partons have been identified 
with current quarks and gluons, i.e. with the fundamental degrees of freedom of the QCD 
Lagrangian. On the contrary, a rigorous derivation of the CQ's from QCD is lacking, but 
CQ's are commonly believed to be quasi-particles emerging from the dressing of valence 
quarks with gluons and quark-antiquark pairs. If CQ's are confined objects, they should be 
connected each other by color strings, which may have their own partonic content. In the 
resolution range in which the sea-quark and gluon content of the strings is not probed, one 
is naturally lead to try to explain the DIS data only in terms of CQ's having a structure. 

The idea to use CQ's as an intermediate step between the current quarks and the hadrons 
is not new at all and indeed it dates back to the seventies [3]. There a two-stage model for 
the parton distributions was proposed, in which any hadron contains a finite number of CQ's 
having a partonic structure. The latter depends only on short-distance (high-Q 2 ) physics, 
which is independent of the particular hadron, while the motion of the CQ's inside the hadron 
reflects the non-perturbative (low-Q 2 ) physics, which depends on the particular hadron. 
Therefore, within such a picture the DIS structure function of a hadron F^(x,Q 2 ) can be 
simply written as the convolution of the structure function of the constituents F$(x / z,Q 2 ) 
with the light-front (LF) momentum distribution fj(z) of the j-th constituent inside the 
hadron H, viz. 



where z is the LF momentum fraction carried by the constituent in the hadron. A convo- 
lution analogous to Eq. (pQ) holds as well for each partonic density in the hadron in terms 
of the corresponding partonic density inside the constituents. The latter can be obtained 
by a deconvolution of available data on a hadron H, provided a reasonable model for the 
wave function describing the motion of the constituents in the hadron H is considered. Then 
the structure function of a different hadron H' can be predicted once its wave function is 
given. Such a procedure has been applied in Ref. [5] to predict the structure function of the 
pion from the known nucleon structure function, and the final result was that the two-stage 
model based on Eq. (JTJ) is supported by data, at least as a first good approximation. 

The following question naturally arises: is the two-stage model a good approximation 
also far from the deep inelastic regime ? In particular, can the model be generalized in such 
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a way to predict hadron structure functions for values of Q 2 below and around the scale 
of chiral symmetry breaking, A x pa 1 GeV ? The aim of this paper is to answer such a 
question by extending the two-stage model in order to include the low-Q 2 regime and to test 
it against recent proton structure function data obtained in Hall B at Jefferson Lab with 
the CLAS spectrometer [Sj. It will be shown that the data exhibit a new type of scaling 
behavior, expected within the generalized two-stage model, and that the resulting scaling 
function can be interpreted as a CQ form factor consistent with the elastic proton (and 
neutron) data. These findings suggest that at low momentum transfer the inclusive proton 
structure function originates mainly from the elastic coupling with extended objects inside 
the proton. We obtain a CQ size of ~ 0.2 0.3 fm. 

The plan of the paper is as follows. The generalization of the original two-stage model 
to low values of Q 2 is presented in Section II and a new type of scaling behavior, which 
should hold for the moments of the structure function, is proposed. In Section III the basic 
theoretical input quantity, i.e. the LF momentum distribution ffiz) of a CQ inside the 
hadron, is discussed and estimated in case of the proton. In Section IV we investigate the 
possible occurrence of the new scaling property in the recent J Lab data 0, as well as the 
possible interpretation of the resulting scaling function as the first experimental evidence of 
the CQ form factor. Our conclusions are summarized in Section V. 

2 Extension of the two-stage model to low momentum 
transfer 

In this Section the original two-stage model of Refs. [H Ej will be generalized in order to 
include the low-Q 2 regime. As a first step, let us develop such a generalization in a simplified 
form, which avoids many complications in the final formulae arising from a complete treat- 
ment of finite-Q 2 effects, but at the same time illustrates the essential physical motivations. 
The proper treatment of kinematical finite-Q 2 effects will be recovered later on in Section 
IV. 

In a DIS experiment at high values of Q 2 the internal structure of a CQ is probed, 
whereas for sufficiently low values of Q 2 such a structure cannot be resolved any more. Gen- 
erally speaking, we expect that the turning point between the high-Q 2 and low-Q 2 regimes 
is around the scale of chiral symmetry breaking, A x pa 1 GeV. As Q 2 decreases below pa A 2 , 
we have two expectations: i) the inelastic coupling of the incoming virtual boson with 
the CQ becomes less and less important, at least because final states are limited by phase 
space effects; ii) the elastic coupling of the incoming virtual boson with the CQ becomes 
more and more important. We point out that at very low values of Q 2 of the order of 
Aq CD [pa 0.1 -f- 0.2 (GeV/c) 2 } the reinteractions among CQ's in the final state, which are not 
considered in our present analysis, cannot be neglected any more (see later on Section 3). 
Therefore, the Q 2 -range where we want to extend the two-stage model is qualitatively given 
by0.1-0.2<Q 2 (GeV7c) 2 < 1-2. 

Let us start by writing the CQ structure function F{ appearing in the convolution formula 



3 



(JTJ) as the sum of two terms F$ = F 2 ^ 4- F^ cl \ corresponding respectively to the inelastic 
and elastic virtual boson coupling with the CQ. Then, the inelastic structure function of a 
hadron, F 2 H (x,Q 2 ), can be written as the sum of two terms, representing the inelastic and 
elastic CQ contributions, respectively. One has 

F 2 H (x, Q 2 ) = £ fdz ff(z) Fi^fi, Q 2 ) +J2 fdz ff(z) if ^(*, Q 2 ) (2) 

J x Z J x Z 

3 3 

where, as previously anticipated, we have kept the simplified convolution form in order to 
avoid up-to-now inessential complications due to finite Q 2 . The elastic part of the CQ 
structure function reads explicitly as 

Ft l) (x',Q 2 ) = G 2 (Q 2 )5(x'-l) (3) 

where 

[G {Q 2 )f = [GUQ 2 )] 2 + t[G> m (Q*)] 2 = [Fi{Q 2 )f + T[ Fi { Q^ (4) 

1 + r 

with Fi(2){Q 2 ) and Ge(m)(Q 2 ) representing the Dirac(Pauli) and electric(magnetic) Sachs 
form factors of the j-th CQ, respectively. Finally, in Eq. ||rE Q 2 /4m 2 with rrij being the 
j-th CQ mass. Thus, the inelastic structure function of the hadron H becomes 

F?(x, Q 2 ) = £ fdz ff(z) F f inel \*, Q 2 ) + E(^(Q 2 )] 2 x ■ ff(x) (5) 

T Jx z 

3 3 

In the DIS regime the elastic CQ contribution is suppressed by the CQ form factors and 
one gets 

F 2 H (x, Q 2 ) ^ DIS £ f dz ff(z) Fl^i*, Q 2 ) (6) 
3 Jx ' z 

On the contrary, for low values of Q 2 the inelastic CQ contribution is expected to become 
negligible and one could have 

F 2 H (x, Q 2 ) -a^< q2 <a, E[^(Q 2 )] 2 x ■ ff{x) (7) 

3 

However, it should be immediately realized that Eq. (j7J) cannot hold at each x value. In- 
deed, at low Q 2 the hadron structure function F^(x, Q 2 ) is characterized by resonance bumps 
emerging over a smooth background, whereas the elastic CQ contribution is expected to have 
a smooth x-shape only, governed by the LF momentum distributions f?{x). Therefore, we 
assume that Eq. (J7j) holds in a dual sense: the x-averages of F^ over each of the resonance 
bumps are representative of the elastic CQ contribution [see the r.h.s. side of Eq. (JTJ)] at the 
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corresponding mean values of x. Such a CQ-hadron duality can be conveniently expressed 
in terms of moments of the hadron structure function, defined as 



M^{Q 2 ) = f 1 dx x n ~ 2 F 2 H {x,Q 2 ) 
Jo 



In a similar way we can define the dual moments as the moments of the elastic CQ contri- 
bution, given by 

M**(Cf) = f Q dx x n ~ 2 £[G,(Q 2 )] 2 x ■ ff{x) (9) 

The occurrence of a CQ-hadron duality for Q 2 < A 2 can be now translated into the domi- 
nance of the dual moments M^ ual (Q 2 ) for low values of n, viz. 

M^{Q 2 )^M^ ual (Q 2 ) (10) 

The limitation to low values of n arises from the fact that as n increases the moment M^{Q 2 ) 
is more and more sensitive to the rapidly varying bumps of the resonances. Therefore Eq. (fTU|) 
cannot hold at very large values of n (see Refs. [7J |HJ El for the case of the parton-hadron 
Bloom-Gilman duality [10J ) . At the same time it should be pointed out that the dual 
relation (fTUj) is expected to hold only for n > 2, because the second moment M 2 (Q 2 ) = 
Jq 1 dx F 2 H (x : Q 2 ) is significantly affected by the low- a; region where the concept of valence 
dominance may become unreliable. 

Let us introduce the squared form factor [F(Q 2 )] 2 defined as 

['Wn = — ^72 — = ^72 ( n ) 

which is normalized to 1 at the photon point. Assuming Si7(2)-symmetric CQ form factors, 
Eq. (JOJ) becomes 

M^ al {Q 2 ) = [F{Q 2 )f-M H n (12) 

with 

M H n = I'dxx"- 1 Yl e]ff{x) (13) 

If one possesses a reasonable model for the CQ momentum distributions fj(x), the moments 
M n can be estimated and therefore the following ratio 

R»(Q 2 ) = M»(Q 2 ) j M H n (14) 

can be constructed starting from the full moments M^(Q 2 ) [Eq. (|SJ)]. The ratio R^(Q 2 ) 
should generally depend on both n and Q 2 as well as on the hadron H. However, when the 
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underlying CQ picture holds, the CQ-hadron duality (fTUj) is expected to hold as well and, 
consequently, the ratio R^{Q 2 ) depends only on Q 2 , i.e. it becomes independent of both the 
order n and the hadron, viz. 



The scaling function, given by the r.h.s. of Eq. (|15|) . is directly the square of the CQ form 
factor, i.e. the form factor of a confined object. The important point is that within our 
generalized two-stage model the new scaling property (|T5|) is expected to occur at low Q 2 . 
We point out that, once the CQ form factor is extracted from known hadron data, the 
moments of the structure function of another hadron can in principle be predicted. 

Let us now introduce the recent results obtained at J Lab [Hj, where the inclusive electron- 
proton cross section has been measured in the nucleon resonance regions for values of Q 2 
below 4.5 (GeV/c) 2 using the CLAS detector. One of the relevant feature of such measure- 
ments is that the CLAS large acceptance has allowed to determine the cross section in a wide 
two-dimensional range of values of Q 2 and x and has made it possible to directly integrate 
all the existing data at fixed Q 2 over the whole significant x-range for the determination 
of the proton moments M^(Q 2 ) with order n > 2. More precisely, the Nachtmann proton 
moments, defined as [TT] 



where r = J 1 + 4M 2 x 2 / Q 2 and £ = 2/(1 + r), have been directly extracted from the data for 
n = 2, 4, 6, 8 As it is well known, the main advantage of the Nachtmann moments (|THj) 
over the Cornwall-Norton moments (JEJ) is that only with the former it is possible to cancel 
out all the finite-Q 2 kinematical corrections due to the non-vanishing mass of the target. 
Thus, in what follows Eq. ([16)1 replaces Eq. (JHJ) for H = p. 

In Fig. ^ the experimental Nachtmann moments M^(Q 2 ), determined in Ref. [B], are 
shown in the Q 2 -range of interest for this work, namely 0.1 < Q 2 < 2 (GeV/c) 2 . The 
contribution arising from the elastic proton peak (x = 1) is not included and therefore, 
from now on, the moments M^(Q 2 ) represent the inelastic part of the proton Nachtmann 
moments. 

The Q 2 -behavior of the moments M^(Q 2 ) shown in Fig.^is characterized by a sharp rise 
at low Q 2 , followed by a smoother behavior for Q 2 > 1 (GeV/c) 2 . However, the dependence 
upon the order n is much more interesting. Indeed, the moments M^(Q 2 ) appear to differ by 
approximately an order of magnitude moving from n to n + 2. As a result, though the range 
of values considered for n is quite restricted (2 < n < 8), the values of the corresponding 
moments are spread over several order of magnitudes. Such a behavior can be qualitatively 
explained within our generalized two-stage model in the following way. Let us assume a very 
simplified and quite rough model for the CQ momentum distribution fj(z) in the proton, 
in which the constituents share exactly just a fraction 1/3 of the proton momentum, viz. 
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Figure 1: Experimental (inelastic) Nachtmann moments M^(Q 2 ) of the proton versus Q 2 from 
Ref- Wi- The dots, squares, diamonds and triangles correspond to n = 2, 4, 6 and 8, respectively. 
The statistical errors are reported, but they are not visible. 



The moments f!13|) simply become 



m: 



n— 1 



(18) 



implying a factor of ~ 1/9 between the orders n and n + 2. Thus, in Fig. |2]we have reported 
the ratio (fT4"j) obtained using the experimental Nachtmann moments M^(Q 2 ) [Eq. (fTEj) ]. 
shown in Fig. CJ and assuming Eq. (fTH|) . It can clearly be seen that with respect to the 
experimental moments M^(Q 2 ) the spread of the ratio R^(Q 2 ) as a function of n has been 
largely reduced. This is an important result obtained with a very simple hypothesis about 
the CQ motion in the proton. Figure 121 shows that there is a clear tendency of the data 
toward a scaling property like Eq. (fTo'j) . 

Any way, we have to consider that Eqs. (|17H18)) imply that the relative motion of the CQ's 
inside the proton is neglected, which is not a reliable assumption in case of light constituents. 
Therefore, in the next Section we perform more realistic estimates of the CQ momentum 
distribution in the proton with the aim of approaching better the scaling property (fTB^) as 
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well as of interpreting the scaling function as a (squared) form factor. 
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Figure 2: Ratio R^(Q 2 ) [Eq. [1$ for H = p] calculated using the experimental Nachtmann mo- 
ments M'P(Q 2 ) [Eq. llty l. shown in Fig. and assuming a delta-like shape for the CQ momentum 
distribution in the proton, namely J2j e< ] fj( x ) = ~ 1/3) [see Eq. \17\) [. The meaning of the 
markers is the same as in Fig. Q 



3 CQ light-front momentum distributions in the pro- 
ton 

Within the two-stage model the basic theoretical input quantity, appearing in Eq. (JT3J), is 
the LF momentum distribution 

= ( 19 ) 

3 

Such a distribution results from the motion of the CQ's inside the particular hadron H and 
in what follows we will explicitly limit ourselves to the case of the proton, which is of interest 
in this work. 
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In order to evaluate the constituent U and D quark distributions in the proton it is 
natural to adopt the Hamiltonian LF formalism ^2]- in terms of the intrinsic LF variables 
£i and ki± (see the Appendix for their definition) the CQ momentum distribution in the 
proton is given by 

= \ £/ [*eAj E - eo <w, |<{^ ; vm}\*?tf m 

where Q = U, D, t v = 1/2, r D = —1/2 and [d^dk^] stands for dki±dk 2 ±dk 3 ± S(ki± + k 2 ± + 
k 31 _) d£id£ 2 d£3 5(£i + £ 2 + £3 - 1)- In Eq. fl23) ^ u p p is the proton LF wave function, whose 
general structure is briefly illustrated in the Appendix, where also all the other relevant 
quantities are defined. Note that the CQ distributions (J2Uj) are normalized as 

[ 1 dzf£ r (z)=2 , f dz &{z) = 1 , (21) 

JO JO 

and satisfy the momentum sum rule 

f dzz[f u (z) + f D (z)] = l- (22) 
jo 

thus, one has 7{z) = [Af^z) + f p D {z))/%. 

In the Appendix the CQ momentum distributions (j2*U|) are explicitly written in terms 
of various SU(6) components characterizing the nucleon wave function [see Eq. ()45|) ]. If a 
completely SU(6) symmetric nucleon wave function is considered, one has always fjj(z) = 
2fp(z) and therefore the LF momentum distribution f P (z) becomes (cf. the Appendix) 

J{z) = [ dk ± dp ± I [dfc] 8(z - 6) ffcft \w s (k,p)\ 2 (23) 

In order to improve the simple delta-like model given by Eq. (|17|) we have calculated Eq. (J23|) 
adopting a gaussian ansatz for the proton wave function ws{k,p), namely 

w s (k,p)oce-^ +3p2 ^ 2 (24) 

where (3 is a parameter. The results of our calculations are reported in Fig.lHlfor various values 
of the CQ mass mu = mo = m keeping the parameter /3 fixed at the value f3 = 0.3 GeV, 
which represents the typical CQ momentum in the proton due to the confinement scale. It 
can be seen that the calculated distribution f P (z) is peak-shaped with a location of the peak 
and a width which sharply depend on m for values of m pertaining to the so-called light-CQ 
sector. The delta-like model (fT7j) . characterized by a zero-width peak located at x = 1/3, 
can be recovered only in the heavy-quark limit m —>■ 00. As the CQ mass decreases, the 
width of the peak increases and the location of the peak moves to values of x less than 1/3. 
Note that: i) the widths are asymmetric around the peaks in order to keep the average 
fraction of the momentum carried by each CQ equal to 1/3 at any value of m, and ii) 
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z 

Figure 3: Light-front momentum distribution fiz) [Eq. \23\) 1. calculated assuming the SU(Q)- 
symmetric gaussian ansatz \2J$ for the proton wave function with (3 = 0.3 GeV. The solid, dashed, 
dotted and dot-dashed lines correspond to a CQ mass equal to m = 0.22,0.33,0.44,0.55 GeV, 
respectively. 

the distribution j^i^z) depends only on the parameter ratio (3/m. Thus, the effects of the 
CQ motion on the shape of f P (z) are very important and should be taken into account, 
particularly for light CQ masses. 

It is well known (see Ref. ^3] and references therein) that a good description of hadronic 
mass spectra requires spin-dependent components in the effective interaction among CQ's. 
Such components generate SU(6) breakings in the proton wave function (see, e.g., Ref. [E]). 
On the contrary the gaussian ansatz (|24j) is a pure SU (6)-symmetric wave function and 
therefore we should investigate SU (6)-breaking effects in the calculation of the CQ light-front 
momentum distribution f P (z). To this end we have considered two of the most sophisticated 
CQ potential models available in the literature, namely the one-gluon-exchange model of 
Ref. ^3] and the chiral model of Ref. ^3] > based on Goldstone-boson-exchange arising from 
the spontaneous breaking of chiral symmetry. The results obtained for f P (z) are shown in 
Fig.|3]and compared with those corresponding to the gaussian ansatz (|24|) for different values 



10 



of the parameter ratio [3/m. It can clearly be seen that, as far as / (z) is concerned, the 
SU(6) breaking contained in the CQ models of Refs. fUJEI] can be approximated to a very 
good extent by using a gaussian ansatz with appropriate values of the parameter ratio (3/m. 




Figure 4: Light-front momentum distribution f (z) [Eq. hl(J\) for H = p], calculated using the full 
proton wave function corresponding to the one-gluon exchange model of Ref. liffl (full dots) and to 
the chiral model of Ref . llftf (open dots). The solid and dashed lines correspond to the case of the 
SU (6) -symmetric gaussian ansatz \21$ with /3/m = 1.2 and 1.8, respectively. 



4 Scaling analysis of the experimental moments 

In this Section we apply our generalized two-stage model to the analysis of the data shown 
in Fig.[H taking into account: i) the motion of the CQ's adopting the gaussian ansatz (|2~H) 
for the proton wave function, as described in the previous Section, and ii) the effects of 
finite Q 2 , which are expected to be relevant due to the Q 2 -range of our analysis [0.1 4- 0.2 < 
Q 2 (GW/c) 2 < 1-2]. 

Let us start by considering the first of the two quoted effects. In Fig. El we have reported 
the results obtained for the ratio Rp(Q 2 ) calculated using the experimental Nachtmann 
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moments M£(Q 2 ) [Eq. (JTHJ)] and assuming the gaussian ansatz (|24j) for the proton wave 
function with (5 = 0.3 GeV and m = 0.25 GeV (corresponding to (3/m = 1.2). The spread 
of the values of the ratio R^(Q 2 ) is drastically reduced with respect to the case of the delta- 
like model (fTTjl [cf. Fig. E]. Note that, as already pointed out (see Section 2) the results at 
n = 2 appear to deviate significantly from those corresponding to larger orders. We have 
checked that the general qualitative shape of the results shown in Fig. El does not change 
significantly when the value of the parameter ratio [3/m is varied. 




I i i i i I i i i i I i i i i I i i i i I 

0.0 0.5 1.0 1.5 2.0 

Q 2 (GeV/c) 2 

Figure 5: Ratio R^(Q 2 ) [Eq. for H = p] calculated using the experimental Nachtmann 
moments MP(Q 2 ) [Eq. Mb)) ] shown in Fig. and assuming the gaussian ansatz \21$ for the proton 
wave function with (3 = 0.3 GeV and m = 0.25 GeV . The meaning of the markers is the same as 
in Fig. Q 

Though the results shown in Fig. El exhibit a drastic improvement toward a significant 
reduction in the dependence of the ratio R^(Q 2 ) upon the order n, the scaling property 
fTS]) is still far from being reached. Moreover, the Q 2 -behavior of R^(Q 2 ) is completely at 
variance with what is naturally expected for a squared form factor. The main drawback is 
clearly the use of Eq. (jlHjl . which is meaningful only at large Q 2 . In our opinion, in order 
to restore a proper behavior of R^(Q 2 ), we have to account for "higher-twist" effects, which 
can be divided into the three following classes: 
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• the inelastic pion threshold, which sets a Q 2 -dependent maximum value for the x-range, 
given by x max = x n = Q 2 /[Q 2 + (M + m n ) 2 — M 2 }. Note that x n largely differs from 1 
at low Q 2 ] 

• kinematical power corrections in the physical region x < x n \ 

• dynamical power corrections due to final state interactions responsible for the resonance 
bumps in the x-space. 

In what follows we will consider the first two effects only. The pion threshold can be simply 
taken into account by multiplying the distribution f P (x) by a threshold factor F t hr(W), 



where W is the produced invariant mass W = y M 2 + Q 2 (l — x)/x, having the property 
Fthr(W < M + m n ) = and F t h r (W — > <x>) = 1. A simple and parameter-free choice 
dictated by pure phase space effects is 



We stress that by means of F t h r (W) we account for that part of higher twists which are 
related to the final-state phase-space constraint. 

The kinematical corrections to Eq. originate from the non-vanishing value of the 
target mass, i.e. the proton mass M. The way to construct such corrections is well known 
in DIS [16j and therefore, for analogy, we replace the distribution f{x) by the quantity 
7tm(£>Q 2 ), S iven explicitly by 

T (PH 2 \ * 2 7(0 6M 2 * 3 rV 



r 3 ^ 2 Q 2 r 4 Jt * 

where ^ is the Nachtmann variable, x = — M 2 ^ 2 /Q 2 ) and ^* = min(l,Q/M) is the 
maximum allowed value of £ (cf. Ref. [9 ). It should be reminded that the value £* is larger 
than the inelastic pion threshold Therefore, the support in which the function 7> M (£, Q 2 ) 
is defined, contains an unphysical region extending from C, n to £*. 

We point out that Eq. expresses the fact that the asymptotic function f p receives a 
series of power corrections having a scale of order of the proton mass M. When the threshold 
factor F t hr{W) is neglected (i.e., F t h r (W) = 1), the use of the Nachtmann moments cancel 
out exactly all the power corrections contained in the r.h.s. of Eq. (|26|). On the contrary, 
when the threshold factor is considered (i.e., F thr (W) ^ 1), only part of the target-mass 
corrections can be reabsorbed by the use of the Nachtmann moments. As a matter of facts, 
for consistency with the experimental data shown in Fig. ^ the Cornwall- Norton moment 
(|T3*|) has to be replaced by a Nachtmann one. In doing that the quantity M n is no more 
independent of Q 2 , and therefore Eqs. (jl2H13j) are now replaced by 

Mt al (Q 2 ) = [F(Q 2 )] 2 -M P n (Q 2 ) (27) 
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with 



Jo 



' n+1 3 + 3(n + l)r + n(n + 2)r 2 
~x~ 3 (n + 2)(n + 3) 



r(l + r) 



(fTM^Q 2 ) F thr (W) 



(28) 



2 



where r(l + r)/2 = dx/d^ arises from the change of variables from x to £. In Eq. (J2%|) we 
have put £* as the upper limit of integration; however, due to the threshold factor (f23J) the 
integration extends only up to £ n and therefore part of the target-mass corrections survives 
after integration. We stress again that this is an important point, because Eq. (J25j) reduces 
exactly to Eq. (fT3*j) when the threshold factor F t h r (W) is disregarded 3 ", in agreement with 
the properties of the Nachtmann moments. 

We have calculated Eq. ()28j) using the target-mass corrected momentum distribution 
/ P (£> Q 2 ) starting from the gaussian ansatz (J2~4~j) for the proton wave function and adopting 
the threshold factor (J2HJ). The results for the moment ratio Rp(Q 2 ) obtained at (5 = 0.3 GeV 
and m = 0.25 GeV are reported in Fig. |U1 It can clearly be seen that the scaling property 
J73J) holds at n > 2 even in a linear scale. Moreover, the scaling function closely resembles 
a squared monopole form factor [F{Q 2 )} 2 = 1/(1 + Tq Q 2 /6) 2 corresponding to a CQ size 
equal to rq = 0.21 fm. 

The quality of the scaling exhibited in Fig. El is extremely good for Q 2 > 0.3 (GeV/c) 2 , 
while it deteriorates at very low values of Q 2 (but still the scaling is approximately satisfied 
within pa 30% even at Q 2 pa 0.1 (GeV/c) 2 ). This finding is not surprising at all and it can 
be understood as follows. Let us consider the Operator Product Expansion (OPE) of the 
moments of the proton structure function in terms of local operators acting on elementary 
(point-like) fields. The so-called higher twists are known to describe correlations among 
partons. Their contribution to the OPE is given by matrix elements of a series of several 
operators O n producing power suppressed terms of the form (A^/ Q 2 )( Tn ~ 2 ^ 2 , where r n is the 
twist and A n is the scale associated with the operators O n . The scale A n is expected to be 
proportional to 1/R n , where R n is the typical average distance of the partonic correlations 
generated by the operators O n . Which kind of higher twists are accounted for by the spatial 
extension of the CQ's ? It is clear that we can distinguish two basic types of partonic 
correlations: those among partons inside the CQ and those between partons belonging to 
different CQ's, which means correlations between CQ's (in the final state). The former 
are characterized by a value of R n close to the CQ size, while the latter correspond to a 
larger value of R n of the order of the confinement (hadronic) size. Correspondingly, the scale 
A n is larger for partonic correlations inside the CQ and smaller for partonic correlations 
among different CQ's. In our model only the first type of higher twists can be thought to 
be accounted for by the CQ form factor in some effective way b . Our model does not include 

a More precisely, when the threshold factor F t hr(W) is disregarded, Eq. reduces at any value of Q 2 
to Jq dx a;" -1 / (x), where x* — min(l, Q/M), as it can be easily checked numerically. 
b Indeed, there is no rigorous derivation of the CQ picture from QCD. 



14 



2 . I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




o I i i i i I i i i i I i i i i I i i i i 

0.0 0.5 1.0 1.5 2.0 

Q 2 (GeV/c) 2 

Figure 6: Ratio K^(Q 2 ) [Eq. for H = p] calculated using the experimental Nachtmann mo- 
ments MP(Q 2 ) [Eq. Mb)) ] shown in Fig. and the theoretical moments M^(Q 2 ) given by Eq. 
The momentum distribution / (£) corresponds to the gaussian ansdtz \2J$ for the proton wave 
function with j3 = 0.3 GeV and m = 0.25 GeV . The dotted line represents the squared monopole 
form factor [F(Q 2 )] 2 = 1/(1 + rq Q 2 /6) 2 corresponding to a CQ size equal to rg = 0.21 fm. The 
meaning of the markers is the same as in Fig. 

power corrections arising from correlations among different CQ's in the final state. Such 
"long-range" higher twists have a low scale of the order of Aqcd, and therefore we expect 
that they should play an important role mainly for Q 2 < Aq CD 0.1 -f- 0.2 (GeV/c) 2 , i.e. 
in the Q 2 -range where the scaling shown in Fig. El is only approximate. The estimate of the 
effects of such "long-range" higher twists is not an easy task and it is well beyond the aim of 
the present paper. Note that the role of the "long-range" higher twists is even more evident 
in the x-space, because these higher twists are responsible for the huge resonance bumps 
which are known to characterize the structure function F$(x, Q 2 ) at low values of Q 2 . 

We should now investigate the impact of different choices of the functional form of the 
threshold factor Fthr{W) as well as of different values of the parameter ratio (3/m. We have 
found that the scaling property ()15j) . clearly exhibited in Fig. EJ is not very sensitive to the 
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specific choice of F t h r {W) and of the parameter ratio /3/m. On the contrary the shape of the 
scaling function is affected both by the choice of F thr {W) and by the value of the parameter 
ratio (3/m. It turns out that: i) the use of the specific form ()25)1 minimizes the scaling 
violation at the lowest Q 2 ; ii) when the ratio f3/m changes from the value 1.2, considered in 
Fig- EE to the value 1.8, the CQ size tq changes correspondingly from 0.21 fm to 0.27 fm. 

We point out that an important consistency requirement can be formulated: the CQ 
form factor extracted from the scaling function and the model used for the wave function 
should be consistent with elastic nucleon data. This is a crucial requirement necessary to 
interpret the scaling function as a (squared) form factor and consequently to get an estimate 
of the CQ size. To check this point we have calculated the nucleon elastic form factors 
adopting the covariant LF approach of Ref . . There the one-body approximation for the 
electromagnetic (e.m.) current operator J M is adopted, viz. 



J* 1 -Ji=J2 



F((Q 2 )r + Fi(Q 2 )^ 

2m 



(29) 



where Q 2 = —q ■ q. The approach of Ref. ^3] is characterized by the choice of a frame 
where q + = 0, which allows to eliminate the contribution of the so-called Z-graph (i.e., 
the pair creation from the vacuum [T7|). The important connection with the Feynmann 
triangle diagram is fully discussed in Ref. JH] and the superiority of the choice q + = for 
the one-body approximation (|2*9~j) is clearly illustrated in Ref. [191 . 
The matrix elements of the (on-shell) nucleon e.m. current read as 



where u(P, vn) is a Dirac spinor, q = P' — P and is the LF wave function of the nucleon 
described in the Appendix, i.e. the same wave function used to calculate the momentum 
distribution / (z). In what follows we adopt a Breit frame where the four-momentum 
transfer q = (q°, q) is given by q° = and q = (q x , q y , q z ) = (Q, 0, 0). 
The nucleon Sachs form factors are then given explicitly by ^3] 

G N E {Q 2 ) = F» (Q 2 ) - ^F?{Q 2 ) = i T r{/+[l - ^a y ]} (31) 



GUQ 2 ) = Fi N (Q 2 ) + F?{Q 2 ) = -^Tr{F^} (32) 

where a y and a z are ordinary 2x2 Pauli matrices. 

We have then calculated Eqs. (|31H32j) using the gaussian ansatz for the nucleon 
wave function and adopting the one-body approximation (|29p with both Dirac and Pauli 
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CQ form factors having the following simple behavior: F((Q 2 ) = ej/(l + Tq Q 2 /Q) and 
Fl{Q 2 ) = + r Q Q 2 /12) 2 . The values of the CQ anomalous magnetic moments, k\j 

and Kd, are fixed by the requirement of reproducing the experimental values of proton 
and neutron magnetic moments. The results of the calculations performed with the same 
parameters adopted in case of the ratio R^(Q 2 ) shown in Fig. El namely (3 = 0.3 GeV, 
m = 0.25 GeV and rq = 0.21 fm, are reported in Fig. [7|as the dashed lines. Note that the 
combination [F(Q 2 )] 2 given by Eq. (JHJ), which is the one relevant for the scaling function 
((EH) , turns out to be almost totally dominated by the contribution of the Dirac form factors 
F((Q 2 ) and it is basically insensitive to the presence of the Pauli form factors F^Q 2 ). 

It can be seen that the calculated form factors slightly overestimate the data, so that we 
can conclude that as a first approximation the scaling function of Fig. El may be interpreted 
as a squared CQ form factor. A better consistency with the data can be reached through 
slight variations of the parameters of our model, namely tq and (3/m. For instance, a nice 
agreement with the elastic data can be simply recovered by increasing the CQ size up to 
rq = 0.33 fm, as shown by the solid lines in Fig. [7| However, we can also ascribe the origin 
of the small discrepancies with the elastic data to the fact that the effects of the dynamical 
correlations among the CQ's in the final state are so far missing in our low-Q 2 model. As 
already pointed out, the inclusion of such effects is not an easy task and it is well beyond 
the aim of the present paper. 

5 Conclusions 

In this work we have first generalized the two-stage model of Refs. [HE], originally developed 
in the DIS regime, to values of Q 2 below the scale of chiral symmetry breaking and above 
the QCD confinement scale, i.e. 0.1 4- 0.2 < Q 2 (GeV/c) 2 < 1 4 2. The essential ingredient 
is the inclusion of the contribution to the inelastic hadronic structure functions arising from 
the elastic coupling at the constituent quark level. We have shown that within such a model 
a new scaling property [see Eq. (fT5|) ] is expected to occur in the inelastic hadronic structure 
functions, provided a reasonable model for the wave function describing the motion of the 
constituents inside the hadron is considered. Moreover, the resulting scaling function can be 
interpreted as the (squared) form factor of the constituent quark, i.e. the form factor of a 
confined object. 

Then we have analyzed the recent experimental determinations of the Nachtmann mo- 
ments of the inelastic structure function of the proton F^x^Q 2 ), obtained at JLab jff, for 
values of Q 2 ranging from ^ 0.1 to w 2 (GeV/c) 2 . The important results we have obtained 
are: 

• the scaling property (fTI)|) is well satisfied by the data; 

• the CQ form factor extracted from the inelastic proton data is overall consistent with 
the one required to explain the elastic nucleon data; 

• the constituent quark size turns out to be ~ 0.2 4 0.3 fm. 
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Our findings clearly suggest that at low momentum transfer the inclusive proton structure 
function F%{x, Q 2 ) originates mainly from the elastic coupling with extended objects inside 
the proton. 

A crucial, mandatory check of the extracted constituent form factor is provided by the 
analysis of the moments of the polarized proton structure function g\(x,Q 2 ). Indeed for 
0.1 -T- 0.2 < Q 2 (GeV/c) 2 < 1 4- 2 a scaling property analogous to Eq. (JT5j) is expected 
to hold also for the Nachtmann moments of gi(x,Q 2 ). The crucial point is that the two 
scaling functions, corresponding to the non-polarized and polarized cases, should coincide 
and provide the same constituent quark form factor. 

Measurements of g\(x, Q 2 ) at low values of Q 2 are still undergoing at J Lab. 
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Appendix: The nucleon light-front wave function 



In this Appendix we briefly recall the basic notations and the relevant structure of the 
nucleon wave function in the Hamiltonian LF formalism (see [EJ). The nucleon LF wave 
function is eigenstate of the non- interacting LF angular momentum operators j 2 and j z , 
where the unit vector z = (0, 0, 1) defines the spin quantization axis. The squared free-mass 
operator is given by 

Ml = ± M ' + m2 (33) 

1=1 si 

where m is the mass of the constituent U and D quarks and 

k± = Pi±-&P± (34) 

are the intrinsic LF variables. The subscript _L indicates the projection perpendicular to 
the spin quantization axis and the plus component of a 4- vector p = (p°,p) is given by 
p + = p° + z ■ p; finally P = (P + , P±) = pi + p2 + P3 is the nucleon LF momentum and pi 
the CQ one. Note that = 1- 

In terms of the longitudinal momentum k iz , related to the variable by 



h z - i f^M - ^[J^ 2 ) (35) 



the free mass operator acquires a familiar form, viz. 

3 3 



M = 5>* = EV™ 2 + N 2 (36) 

i=l i=l 

with the three-vectors fcj defined as 

ki = {ki±, kiz) (37) 

Note that ki are internal variables satisfying ki + k 2 + k 3 = 0. Disregarding the color variables, 
the nucleon LF wave function reads as 



<{6^;i4t,W> = J§§YfJ: ■ ({k;vm}\x7) (38) 

V 0slS2s3 {jyj}. 

where Vn is the third component of the nucleon spin, the curly braces { } mean a list 
of indexes corresponding to i — 1, 2, 3, and {ji) is the third component of the CQ spin 
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(isospin). The rotation 7tf , appearing in Eq. (|S5|l. is the product of individual (generalized) 
Melosh rotations, viz. 



3=1 



J- 



m) 



(39) 



where 



Rj(kj±,£j,m) 



m + £ 3 -Mo — ia^ ■ (n x fejjj 
m + ^M ) 2 + |fc i± | 2 



(40) 



with c? being the ordinary Pauli spin matrices. 

Neglecting the very small P- and D-waves in the nucleon (cf. |14j ) we can limit our- 
selves to canonical (or equal-time) wave function corresponding to a total orbital angular 
momentum equal to L = 0; one has 



1 

{{h;ViTi}\xN) = w s(k,p)-j= 



ai 



$ 01 + $ 



ao 



+ w A (k,p) 



V2 



01 

V N T N 



10 

V N T N 



(41) 



where Ws(k,p), ws> B (k,p), ws> a (k,p) and w A (k,p) are the completely symmetric (5), the 
two mixed-symmetry (S 1 ^ and S^) and the completely antisymmetric (A) wave functions, 
respectively. In Eq. (j4~Tj) the variables fc and p are the Jacobian internal coordinates, defined 

as 



fe 



V 



h - k 2 



2k 3 - (k t + fee 



(42) 



with fej given by Eq. (J37|) . Finally, the spin-isospin function ^j!^ 2 ) corresponding to a total 
spin (1/2) and total isospin (1/2), is defined as 



$ ™ = E<^4^i 2 M s > (£ 12 Ms^ 3 |~^> 



"2 

1 

2 



(43) 
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where Su (T12) is the total spin (isospin) of the quark pair (1,2). The normalization of the 
various partial waves in Eq. (|41|) is 



dkdp Ws(k,p) 
dkdp ws' s {k,p) 
dkdp WA(k,p) 



dkdp w S ' a (k,p) = Ps>/2 



Pa 



(44) 



with P s + P s > +P A = 1. 

Disregarding the completely antisymmetric component WA{k,p), which is usually quite 
negligible in the nucleon (cf. [O]), the constituent U and D momentum distributions, defined 
in Eq. (j2~U]l . read explicitly as 



dkj_dp± J [d£i] 5(2 - £1 



\ws(k,p)\ 2 + \ws' s (k,p}\' 



dk ± dpj_ / [d£i] 5(z — £i 



+ 



EiE 2 E 3 

\ w S' a (k,P)\ 2 + w s (k,p) w S ' s (k,p) 
E\E 2 E 3 



\ws(k,p)\ 2 + \w s >(k,p)\ 2 + \w S ' a (k,p)\ 2 - 2w s (k,p) w S ' s (k,p) 



(45) 



It can be seen that the relativistic composition of the CQ spins (i.e. the Melosh rotations) do 
not affect at all the (unpolarized) LF momentum distribution fn(z). Moreover, any flavor 
dependence of fn(z) turns out to be driven by the interference between the completely 
symmetric (S) and mixed-symmetry (S' s ) wave functions, the latter being generated mainly 
by the spin-spin component of the interaction among CQ's which are present both in the 
one-gluon-exchange model of Ref. jTHj and in the chiral model of Ref. [T^j. In the limit of 
exact SU(6) symmetry one has ws' = ws> a = and Eq. is recovered. 
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Figure 7: Elastic Sachs form factors of the nucleon, G P E (Q 2 ) (a), G P M (Q 2 ) (b), G%(Q 2 ) (c), 
—G% I {Q 2 ) (d), calculated using the covariant LF approach of Ref. \14\l - The gaussian ansatz \2J$ is 
adopted for the proton wave function with (3 = 0.3 GeV and m = 0.25 GeV. In the one-body current 
both Dirac and Pauli CQ form factors are included, namely: F((Q 2 ) = ej/(l + Q 2 /6) and 
F^iQ 2 ) = + Tq Q 2 /12) 2 . The dashed and solid lines correspond to rQ = 0.21 and 0.33 fm, 

respectively (see text). The values of the CQ anomalous magnetic moments, kjj = —0.064 and 
K£, = 0.017, have been fixed by the requirement of reproducing the experimental values of proton 
and neutron magnetic moments. In (a) and (b) full and open dots, open and full squares are 
the experimental data from Ref. \2(tij(a-d). respectively. In (c) open squares, full squares, open 
diamonds, open triangles, full dots, full diamonds and triangles are the data from Ref. [21}(a-g), 
respectively. In (d) full dots, open dots, full and open squares are the data from Ref. ' L 22](a-d), 
respectively. 
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